Black Holes with Skyrme Hair 



Noriko Shiiki and Nobuyuki Sawado 
Department of Physics, Tokyo University of Science, Noda, Chiba 278-8510, Japan 

February 7, 2008 



Abstract 

This paper is intended to give a review of the recent developments on black 
holes with Skyrme hair. 

The Skyrme model is an effective meson theory where the baryons are iden- 
tified with topological solitons, so-called skyrmions. The baryon number B 
corresponds to the topological charge. Thus the model gives a unified descrip- 
tion of hadron physics in terms of meson fields. The spherically symmetric 
B = 1 skyrmion and axially symmetric B — 2 skyrmion were found numeri- 
cally. Upon collective quantization, they produce correct nucleon and deuteron 
observables respectively. Approximate solutions with higher baryon numbers 
were also constructed by using the rational map ansatz. Depending on the 
baryon number, they exhibit various discrete symmetries with striking similar- 
ities to BPS monopoles. 

The Einstein-Skyrme system is also known to possess black hole solutions 
with Skyrme hair. The spherically symmetric black hole skyrmion with B = 1 
was the first discovered counter example of the no-hair conjecture for black 
holes. Recently we found the B = 2 axially symmetric black hole skyrmion. In 
this system, the black hole at the center of the skyrmion absorbs the baryon 
number partially, leaving fractional charge outside the horizon. Therefore the 
baryon number is no longer conserved. We examine the B = 1,2 black hole 
solutions in detail in this paper. 

The model has a natural extension to the gauged version which can describe 
monopole black hole skyrmions. Callan and Witten discussed the monopole 
catalysis of proton decay within the Skyrme model. We apply the idea to the 
Einstein-Maxwell-Skyrme system and obtain monopole black hole skyrmions. 
Remarkably there exist multi-black hole skyrmion solutions in which the grav- 
itational, electromagnetic, and strong forces between the monopoles are all in 
balance. The solutions turn out to be stable under spherically symmetric linear 
perturbations. 
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1 Introduction 

Application of the Skyrme model to the physics of the early Universe or equivalent 
high-energy physics is an interesting subject. 

The attempts to explain (anti-)baryon production in high-energy collisions within 
the Skyrme model were made in Refs. 0I2HB]- The scenario for the production is 
the same as the one for producing other topological defects such as monopoles, cos- 
mic strings or domain walls in the early Universe via the Kibble mechanism [U IS] . 
Those defects are formed after a phase transition in which global symmetry is bro- 
ken spontaneously. In the case of skyrmions, broken is the chiral symmetry and 
as a result one of the possible orientations of the chiral field is chosen randomly in 
the internal symmetry space to form the defects. The use of the Skyrme model is 
expected to describe nonperturbative phenomena to which perturbative QCD is not 
accessible. 

Another interesting application of the Skyrme model is the monopole catalysis 
of proton decay. The monopole catalysis of proton decay was studied firstly at the 
lepton and quark level by Callan jSJ, Rubakov [Zj and Wilczek [B] independently. 
In the grand unified theory there was a symmetry between baryons and leptons at 
energy scale 10 15 GeV and therefore it would have been possible for the following 
reaction to be driven through the anomaly 

p + Monopole — > e + + pions + Monopole . (1) 

Interestingly this reaction occurs at a strong interaction rate unsuppressed by any 
small-coupling-constant effects. Callan and Witten proposed not quark-monopole 
but proton-monopole interactions by using the Skyrme model ,9 to make a more 
realistic estimation of the catalysis cross section possible. The decay process can be 
explained because in the presence of a monopole the baryon number is no longer 
equal to the topological charge of the meson field. In fact, there exist non-topological 
solitons with non-zero baryon number which can decay without topological prob- 
lems. The extension to non-abelian monopole catalysis was also studied in Ref. |lflj . 

For the study of the high-energy physics in the early Universe, it may be im- 
portant to consider the effects of gravity on baryons. With the present value of the 
gravitational constant, those effects are insignificant upto the Planck energy 10 19 
GeV. In fact, in the Einstein-Skyrme theory, the Planck mass is related to the pion 
decay constant F n and the coupling constant a by M p i = F^^/An/a. To realize the 
realistic value of the Planck mass, the coupling constant should be extremely small 
with a ~ O(10~ 39 ), which makes the theory little different from the theory without 
gravity. However, some theories such as scalar-tensor gravity theory and Kaluza- 
Klein theory predict the time variation of the gravitational constant I12[ 113] . 
And also theories with extra dimensions predict that a true Planck scale is of order 
a TeV. Thus there may have been an epoch in the early Universe where the gravi- 
tational effects on baryons were significant. We consider those effects worth being 
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studied in the Skyrme model. The advantage of the Skyrme model to more realistic 
nucleon models is that it is straightforward to couple to gravity. Although inves- 
tigating the Einstein-Skyrme system is not expected to provide any quantitatively 
reliable results, we cannot exclude the possibility that it may give some qualitative 
insight into the gravitational interaction of baryons. 

The Einstein-Skyrme system has been studied by various authors. The first ob- 
tained solutions in this system are spherically symmetric black holes with Skyrme 
hair |141 1151 HE] . Later, regular solutions for B = 1 |151 1161 IT7| and axially sym- 
metric black hole and regular solutions for B = 2 JH] were found. The extended 
models to SU(3) and SU(N) were also studied in Refs. jTS]. It has been observed 
that microscopic black holes can support Skyrme hair admitting fractional baryon 
number outside the horizon. This configuration can be interpreted as a skyrmion 
partially absorbed by the microscopic black hole. Therefore the model provides a 
semiclassical framework to study the absorption rate of a proton by a black hole of 
comparable size. However this process is rather insignificant because black holes of 
the size of a proton have large fluxes of Hawking radiation |20] • Situations in which 
baryon decay process might become more realistic and significant over the Hawking 
radiation occurs when the black hole carries electric or/and magnetic charge with 
which the skyrmion interacts electromagnetically as well as gravitationally. Since 
a charged black hole in general has an effective temperature lower than that of a 
Schwarzschild black hole of the same mass, the Hawking radiation effect is less. 
Especially interesting is the extremal black hole which has a vanishing effective 
temperature, so the Hawking radiation may even vanish |21| . According to this 
speculation, we analyzed monopole black hole solutions with Skyrme hair j^j. This 
model provides the semiclassical framework in which to study monopole black hole 
catalysis of proton decay. Although macroscopic charged black holes perhaps do 
not exist in nature, microscopic charged black holes may have been created in the 
early Universe and remain as stable relics today. Indeed the GUT monopole has a 
typical radius Rm ~ MZ and a mass g~ 2 Mx, the Schwarzschild radius for the mass 
is then Rbh ~ 2Gg~ 2 Mx and thus Rm/Rbh ~ 2Gg~ 2 M x where g 2 is the cou- 
pling constant renormalized at the X-boson mass scale Mx- Hence for a sufficiently 
large magnetic charge p = 1/g and heavy gauge boson mass, the monopoles could 
have undergone gravitational collapse to form monopole black holes. This monopole 
black hole has an additional internal degree of freedom for electric charges being a 
dyon black hole analogous to the 't Hooft-Polyakov monopoles being the Julia-Zee 
dyon [2*H] . 

In subsequent sections we review the recent developments on black holes with 
Skyrme hair. Section 2 contains an introduction of the Skyrme model along with 
its brief historical review. In section 3 we review the B = 1 black hole skyrmion 
solution and its stability analysis. In section 4 the Einstein-Skyrme model coupled 
to abelian gauge fields is studied to obtain B = 1 monopole black hole skyrmion 
solutions. Its stability is examined in detail and the technique to study monopole 
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black hole catalysis of baryon decay is discussed briefly B = 2 black hole skyrmions 
are examined in section 5. It is shown that the energy density and baryon density 
of the solutions are torus in shape having a horizon at the center. Conclusion and 
Discussion are in section 6. 

Throughout this paper we use the metric signature (—,+,+,•••) and the Ein- 
stein summation convention. Greek indices are used to denote spacetime compo- 
nents of a tensor, while Italian indices are used to denote purely spatial components. 
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2 The Skyrme Model and Skyrmions 

In this section we give an introduction of the Skyrme model, reviewing its historical 
development. A more detailed and complete review of the Skyrme model is provided 
by Zahed and Brown in Ref. |24| . 

2.1 Historical Review 

It has been known that in the large-7V c effective meson theory, baryons emerge as 
solitons. This recognition has a long history. 

In the Standard Model QCD is the theory for strong interactions. However it 
is SU(3) gauge theory and there is no hope of solving it exactly and therefore an 
approximation method is necessary for the qualitative estimation of QCD. The most 
convenient approximation method would be a perturbation which, in fact, worked 
very well for QED. However it seemed that there was no small parameter to expand 
in low-energy QCD since the coupling constant is no longer small at low energies. 
The unexpected possible parameter was proposed by 't Hooft [2^]. Although in 
the real world hadrons have definite color degrees of freedom N c = 3, he suggested 
considering the number of colors N c as a free parameter and obtaining the series 
expansion in powers of 1/N C . Then a detailed analysis of the diagrams shows that 
the limit of N c — > oo fixing g%N c is equivalent to the effective chiral field theory for 
mesons where 53 is a coupling constant. Thus 't Hooft succeeded to prove that the 
chiral effective Lagrangian is indeed low energy QCD. The next question was, how 
can baryons be incorporated in this large N c effective meson theory? 

Witten gave a rigorous answer to this question 26 . In his paper he showed that 
baryon mass is of order N c which is inverse of the expansion parameter 1/N C . Hence 
in the large- N c limit, they become singular. This is a typical feature seen in solitons. 
From this fact, Witten conjectured that baryons are solitons in the large- N c meson 
theory. 

Much earlier than Witten, in the flow of development of the effective theory of 
QCD, Skyrme attempted to construct a unified theory of strong interactions based 
on the meson field alone [27]- He explicitly constructed the minimal Lagrangian if the 
non-linear sigma model in d = 1 + 3 that supports topological soliton solutions and 
suggested that they are baryons |28| . However this remarkable idea was neglected 
until the 80s. 

A unification of the ideas came in 1983. Witten studied solitons in the current 
algebra effective Lagrangian, which is the SU(N C ) chiral models with the QCD 
anomalous term (Wess-Zumino term [301. In QCD, the decay process KK — > 
7T + +7r~+7r° is anomalous in the sense that it is observed in nature but is forbidden in 
the chiral model because of the symmetry U — * which conserves the parity of the 
number of mesons. By including the Wess-Zumino term in the action, the anomalous 
process can be correctly described [H^. Witten showed that the anomalous piece 
of the current is responsible for the soliton solutions and their quantum numbers 
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such as spin and baryon number. In particular for the N c = 2 case which is the 
chiral model with N c flavor degrees of freedom, the Wess-Zumino term vanishes. 
However from the homotopy argument 7^(577(2)) = Z, there are still solitons. 
Besides Tt4,{SU{2)) = Z2 suggests that in a suitably compactified space-time, there 
are two topological classes of maps from space-time to SU(2) and hence the solitons 
can be spin-^ particles jSlj. A detailed analysis for the property of the skyrmion as 
a nucleon was performed in Ref. [33] u P° n quantization of its collective coordinates. 



2.2 The Model 

In the lowest order, the chiral Lagrangian is given by 

C = ^ tr (U^d^UU^d^U) , U = e 2i ^ F * (2) 

in the exponential parameterization where F n = 186 MeV is the pion decay constant, 
t are the Pauli matrices and tt are pion fields. The static energy is then given by 

f F 2 

E= d 3 x^tv (U^diUU^diU) . (3) 
J 16 

One can show that this model does not support a topological soliton as follows. Let 
us introduce the dimensionless variable x defined by x = ax. Then for a static 
configuration, the energy can be written as 

/F 2 
d 3 x-^ti(U^d t UU^d t U). (4) 

The integrand on the right-hand side is non-negative and hence the energy is mini- 
mized at a = with E = 0. Therefore the solution is trivial. 

Skyrme introduced a new term in the fourth order derivative which retains the 
chiral symmetry of the model but supports a soliton solution [2B]- The so-called 
Skyrme model is defined by 

F 2 1 
Cs = if tr ^ 9 ^ Uu1 W + ^2 tr \% UU \ duUrf} 2 (5) 

where a is a dimensionless parameter whose value can be fixed by experiment. The 
static energy is then given by 

e = J d^^tviu^diUu^di^ + ^ti^uu^djUU^Y (6) 

One can see that this model indeed supports a nontrivial soliton solution by following 
the earlier argument. Expressing the static energy in terms of the rescaled variable 
x, one obtains 

E = a [ d 3 x^tr (U^diUWdiU) + - [ d 3 x — tr \diUU\ djUU*] 2 . (7) 
J 16 a J 32a 2 
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Since both integrands on the right-hand side are nonnegative, the energy is mini- 
mized at 

Let us find a topological soliton in the Skyrme model. We introduce = U^d^U 
and write the energy in terms of An 



E = J d 3 xtv l^AiA] + ^^(e ijk AjA k ){eii m AiA m yX . 



(8) 



The boundary condition for E to be finite is 

Ai — ► as |x| — > oo (9) 

which is equivalent to saying that U approaches some constant matrix at infinity. 
Without loss of generality, we define this constant as the unit matrix 

U^I as|x|^oo. (10) 

Then the spacetime is compactified to the three-sphere S 3 . 

There is a lower bound for the energy (JSj). From the Cauchy-Schwarz inequality, 

K „ 1 x 2 



one can obtain 



-^A i -—e ijk A j A k ) >0, (11) 



E>^-jd 3 x |tr (eijkAiAjAk)] . (12) 
The topological current of the skyrmion is 

f fiupcr 

B» = -^tr(A„A p A a ) (13) 

with d^B^ = 0. The topological charge is given by the zeroth component of the 
current 

B = Jd 3 xB° = --L J d^trieijkAA.Ak) (14) 

which corresponds to the winding number of the map S 3 — ► S 3 and characterized 
by an integer as tt^^S 3 ) = Z. In the Skyrme model, topological charge is identified 
with the baryon number and hence skyrmions are interpreted as baryons |28j . 
From Eqs. (|12|) and 1)14(1 one finds the Bogomol'nyi bound 

E>^L\B\. (15) 
a 

At present, no soliton solution saturating this bound is found. But the skyrmions 
which have been obtained numerically probably represent the global minimum of 
the energy for given B. 
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Figure 1: Profile / as a function of radial distance x = aF n r. 

The Skyrmion can be found by using the hedgehog ansatz 

U = e i/(r) *- r = cos f{r) + in ■ f sin f(r) (16) 
where n = x/r with the boundary conditions 

/(0)=7T, /(oo) = 0. (17) 
Inserting the ansatz into the energy functional (jBJ), one gets 

„2 



E = 4tt drr 



o 



f/'2 + 2sin2 / ^) + 1^ f sin 2 / + 2//2 
8 \ r 2 ) 2a 2 r 2 \ r 2 



(18) 



where a prime denotes derivative with respect to r. The static solution can be given 
as an extremum of the energy. Hence it satisfies 5E/5F = 0. Equivalently 

%2 , o • 2 ,\ -// . x f i . w2 • o - sin 2 /sin 2/ 
— + 2 sin / 1 / + -/ + / sin 2/ = (19) 

where we introduced a dimensionless variable x = aF n r. The solution of the equation 
(|19l) satisfying the boundary conditions (|17|) is shown in Fig. ^ For this solution, 
we have 

B = -^[2/-sin2/]°=l. (20) 

Thus this solution represents a nucleon. For higher baryon numbers, B = 2 skyrmions 
were obtained numerically and shown to be axially symmetric |34|I35|. Braaten et al. 
constructed skyrmions upto B = 6 by descretizing the model on a cubic lattice 36 . 
Interestingly, it has been shown that multi-skyrmions with B > 2 exhibit various 
discrete symmetries analogously to multi-BPS monopoles in the use the rational 
map ansatz proposed in Ref. |37j . 
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3 B = 1 Black Hole Skyrmions 

As stated earlier, it is straightforward to couple the Skyrme model to gravity. The 
Skyrme Lagrangian and topological current are written in a covariant manner. After 
imposing the suitable ansatz and boundary conditions on both the chiral field and 
metric, one can solve numerically the Einstein equations coupled to the chiral field 
to obtain black hole solutions with Skyrme hair. The Einstein-Skyrme system was 
firstly studied by Luckock and Moss Jl] where the Schwarzschild black hole with 
Skyrme hair was obtained numerically. This is a counter example of the no-hair 
conjecture for black holes [HE]- They observed that the presence of the horizon 
in the core of skyrmion unwinds the skyrmion, leaving fractional baryonic charge 
outside the horizon. The full Einstein-Skyrme system was solved later to obtain 
spherically symmetric black holes with Skyrme hair |15U16j and regular gravitating 
skyrmions |16l I17j . In this section we review the B = 1 black hole skyrmion solution 
and analyze its stability following Refs. fPHTnl [TCI IT71 133] . 

3.1 Field Equations and Static Solutions 

The Einstein-Skyrme system is defined by the Lagrangian 

C = C G + Cs (21) 
li + §<Ttr {L,L V ) + A-gTsTU L p ] [L u , L CT ]) 



16vrG 16 a v ^ 32a 2 " 
For B = 1, let us impose the hedgehog ansatz on the chiral field 

U(r) = cos f(r) + in ■ f sin f(r) . (22) 

Correspondingly, we shall impose the spherically symmetric ansatz on the metric 

ds 2 = -N 2 (r)C(r) dt 2 + -J--dr 2 + r 2 dtt 2 (23) 

C(r) 

where we have defined 

2Gm(r) 



C{r) = 1 

At the horizon r = rh, we have C(rh) = 0, that is, m{rh) = rh/(2G). Inserting 
these ansatz into the Lagrangian (|2*T|) . one obtains the static energy density for the 
chiral field 

where we have introduced dimensionless variables 

x = aF n r = r/0.28fm , rh{x) = aF n Gm{r) . 
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It should be noted that in the presence of a black hole, the matter field is defined 
only outside the horizon and therefore the integral over the space is performed from 
the horizon to infinity. 

The covariant topological current is defined by 



=tr (U^dvUU^dpUU^daU) . 



24vr 2 yj-g 

whose zeroth component corresponds to the baryon number density 

1 1 /'sin 2 / 



2vr 2 N r 2 



We impose the boundary conditions on the profile function 

/(oo) = 0, 



(25) 



(26) 



(27) 



which determines the value at the horizon /(r/j) = fh- Then the baryon number 
becomes 



B= [ ^B°d 3 x = -- /°sin 2 M = ^(2A-sin2/ h ). 



(28) 



This shows that the solution possesses fractional baryonic charge when fh < it. 

The field equations for the gravitational fields N(x) and m(x) can be derived 
from the Einstein equations 



which read 



N' 



■m 



a 

4 1 ' + 



8 sin 2 F 



Nf 



12 



a 



(x 2 + 8 sin 2 f)Cf' 2 + 2 sin 2 / + 



(29) 

(30) 
(31) 



where we have defined the coupling constant a = 4nGF 2 . The variation of the static 
energy (|24[) with respect to the profile f(x) leads to the field equation for matter 



NC (x 2 + 8sin 2 f) 



-(x 2 + 8sm 2 f)N'Cf + 1 + 



4 sin 2 / 



+ 4Cf 



r2 



xiVsin2/-2(x + 4sin2// / )iVC/ / - 2 ( 1 + 



8 sin 2 / 



(fh — m'x)Nf' 



■ (32) 



To solve these coupled field equations, let us consider the boundary conditions for 
the gravitational fields. Expanding the fields f(x), rh(x), N(x) around the horizon 
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x = Xh and substituting into the field equations ()3(J|) - (|32|) . one obtains upto second 
order 

f{x) = fh + fi(x-x h )+0((x-x h ) 2 ) 
m(x) = + rhi(x - x h ) + 0((x - x h ) 2 ) 
N(x) = Nh + N^x-x^ + Oiix-Xhf) 

where fh and Nh are shooting parameters which should be determined so as to 
satisfy the boundary conditions at infinity /(oo) = and iV(oo) = 1, and 

a ( . 2 2sm 4 /fc\ 

4 [ Sm fh + ) (33) 

{x\ + 4 sin 2 f h ) sin 2 f h 
x h (x 2 h + 8sm 2 f h )(l-2m 1 ) 
a ( 8 sin 2 f h \ 2 

i + ) Nhfi • (35) 

The dependence of the profile function on the coupling constant and the horizon size 
is shown in Fig. |2j There are two branches of solutions for each value of the coupling 
constant. Let us call the solution with larger (smaller) values of fh as upper (lower)- 
branch. The skyrmion shrinks as a or/and Xh increase for upper-branch while it 
expands for lower-branch. Fig. shows the dependence of the value of the profile 
function at the horizon on the horizon size. One can see that for a^O, all black 
hole solutions converge to globally regular solutions as n ^ 0. As a approaches 
to zero, the lower-branch solutions converge to the Schwartzschild black hole with 
Skyrme hair. On the other hand, the upper-branch solutions converge to the n = 1 
colored black hole solution [10] . The maximum value of the coupling constant with 
which black hole solutions exit is 0.126. 



rhi = 

fl = 
Nt = 
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Figure 2: Profile function f(x) with a = 0.00,0.10 and Xh = 0.05,0.10 for upper 
and lower branch. 




Figure 3: The shooting parameter fh as a function of the horizon size Xh with 
a = 0.00, 0.05, 0.10. Solid line shows the stable branch and dotted line the unstable 
branch. 
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3.2 Linear Stability Analysis 

The linear stability of the B = 1 black hole skyrmion was studied in detail in 
Refs. [391 117j . In order to examine the stability, let us consider the time-dependent 
Skyrme action given by 



vre 2 Fi 



(- 



1 



f + Cf ,2 )u + v 



e d C 



e s dx 



where we have defined 

5 = logiV , u = x 2 + 8 sin 2 / , v = sin 2 fix 2 + 2 sin 2 /) 



(36) 



(37) 



Then time-dependent field equation can be obtained by taking variation with respect 
to / as 



(e«Cuf'y + - 



1 



2 \e s C 



f 2 ~ e s Cf' 2 ) u f 



e 5 v f 



1 



uf 



(38) 



where Uf = Su/5f and Vf = 5v/5f. From the action (|36|) . the time-dependent 
Einstein equations are derived as 



Goo = 87rGToo 
Gn = 8ttGTh 



l- C -C'x = ^- 
4 



' f 2 + cAu + % 



(e 2S C )' 

P 25 



-1 + C + 



e 2S C 



of 



-x 



1 



P + Cf 2 )u-^ 



,26 C 



&*> 4 

which read to the following two equations for the gravitational fields 

J2 / 1 



2i< 

2 



X 



Ax \e 2S C 2 



f 2 + f >2 )u 



-(Cx)' + 1 = ^v + C5'x. 
2x z 



(39) 
(40) 

(41) 
(42) 



Let us consider the small radial fluctuations of the profile and gravitational fields 
around the static classical solutions denoted by /o, Nq and Cq as 



f(x,t) = Mx) + fi(x,t) 
5(x,t) = Sq(x) + 5i(x,t) 
C(x,t) = C (x) + Ci(x,t). 

Substituting into Eqs. ()41j) and (|4*2*)l gives the linearized equations 



(43) 
(44) 
(45) 



5' 1 = ^-{2uM+u f Xh) 

2 

-(e 6 °Cix)' = ^e^vfji + e 5o C 5'iX . 



(46) 
(47) 
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Eq. (|46|) and the classical field equation derived from Eq. (|38j) 



^ = (e s °C u f>y - U 5 °c u f x 

are inserted into Eq. (|47|) and resultantly one gets 

-(e s °C lX )' = ^(e 5 °C u f>f i y 



which can be integrated immediately 



2~r 



Co^o/o/i • 



(48) 



(49) 



(50) 



Similarly we shall linearize the field equation (|38|) . Using Eqs. (|4fijl. (jl8|) and (|50|) 

one arrives at 

1 



fi 



(51) 



where 



-( e 5 oc oU/o ^)'+ (^e s °C Q ulf> 2 ^ -^C u u f X 



+ 2 e G o«//o/o + x 2 • 
Setting /i = ^(x)e ia; V\/^o , Eq. dHU) becomes 



(52) 



(e*°CoO' + 



1 



+ — E7b 

■U 



* 2 1 t 



Let us introduce the tortoise coordinate x* such that 

dx* 1 
dx ~ e*C 

with — oo < x* < +oo. Eq. (j53[) is then reduced to the Strum-Liouville equation 

d\ 



(53) 



(54) 



*2 



where 



Uo = e do C 



e So C n U ° 



+ —u 

u 



(55) 



(56) 



If the black hole skyrmion is stable, this equation has no negative mode which 
induces exponential grow in £. Numerically solving the Eq. (|53jl or (|56j) . one finds 
no negative mode in the upper branch and one negative mode in the lower branch. 
Therefore, it is concluded that the black hole skyrmion in the upper branch is stable 
and unstable in the lower branch. 
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4 Monopole Black Hole Skyrmions 

In this section, charged black hole solutions with Skyrme hair are discussed. We find 
stable non-topological skyrmion solutions as well as topological ones in a background 
of charged black holes. The Skyrme model is valid at a greater scale than the quark 
confinement scale where the gauge fields can be assumed to be abelian. Although 
more complicated non-abelian monopole solutions are possible, depending on the 
details of the Higgs sector we shall restrict attention to the simplest case. The 
black hole mass for abelian monopoles has the lower bound 

M = -j= = pM pl « 2.54 x 1(T 7 Kg , (57) 

where we denoted the magnetic charge p ~ 11.7 and the Planck mass M p i ~ 2.1768 x 
10~ 8 Kg. We shall also examine the stability of our solutions. The main obstacle 
to proton decay around the monopole black hole is electric charge conservation. 
The black hole cannot swallow the proton whole because this would tip it over the 
extremal limit. It should come as no surprise, therefore, that the non-topological 
solutions are stable. The introduction of charged fermions, which can carry the 
electric charge away, will be also briefly discussed. 

This work was carried out with the collaboration of I. Moss and E. Winstanley. 



4.1 Field Equations and Static Solutions 

The Lagrangian is based on a gauged version of the original Skyrme Lagrangian 
constructed by Callan and Witten ,9 a . The natural extension to the charged SU{2) 
chiral field is gauging the Skyrme model. For a gauge transformation 

U -► U + iea [Q, U] , A fl ^A fl + d^a , (58) 

one can define the covariant derivative 

DpU = d^U - ieA^ [Q, U] (59) 

where A^ is the photon field, e is the charge of proton electric charge in unrational- 
ized units, and 

Q=(| (60) 

is the charge matrix of quarks. Replacing merely the derivatives in the Lagrangian 
with the covariant derivatives DJJ = d^U — ieA^ [Q, U] is, however, not sufficient in 
the sense that QCD anomalous processes, such as 7r° — > 77 or the 77r + 7r~7r° vertex, 
are not included. This is a manifestation of the chiral symmetry breaking due to 
the presence of the U(l) gauge fields. Thus the correct Lagrangian should include 
an anomalous term Ca- The anomalous term arises from the Wess-Zumino term 
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which vanishes in the SU(2) Skyrme model. It has been worked out in Refs. [HI 131 j 
that the gauge invariant form of the Wess-Zumino term T(U) has additional terms 
given by 

f (U, A M ) = T(U) + ejd 4 x A p J» - ie 2 J d 4 x F pu A p T a (61) 

which reads 

C A = T^2^ Upa Mm tr {QLvL p L a - QR v R p R a ) - ie 2 F pv A p T a ] (62) 

where the left current L p = U^d^U and right current = Ud p U^ and 

T a = tr (Q 2 L a - Q 2 R a + ^QU j QUL a - -QUQU^R a ) (63) 

have been defined. Consequently, we have the following gauged Einstein-Skyrme 
Lagrangian 

C = C S + C A + Cem + C G (64) 

where 

Cs = ^tv(U^D p UU^D p U) + ^ti([U^D p U,U^D u U} 2 ) 

c A = -^e^p°A^\Q{u^d v uu^d p uu^d a u + d v utfd p utfd a utf)} 

IP 

x tr [Q 2 (d a U)tf + Q 2 U\8 a U) + ~Q(d a U)Qtf - ^QUQ(d a U^)} 



and the free actions are 

£em = ~16^ F ^' Cg= 16^G R - 



The gauge-invariant baryon current is then given by 



e 



pupa 



Bfl = ^^\^(LvLpLo)+ZieA v X,YQ(L p L u -R p R u ) + ?>ied v A p tYQ(L a -R u ) 



24tt 2 

c pvp<T e pvpa 



} j _ 2 tr (L v L p L a ) + — - d u [ZieApti Q(L a - R a )} (65) 

which shows the additional current is a total divergence. 

In the spherically symmetric case with a magnetic charge the gauge field has the 
form 

A = p(l - cos 6) dcf> + <5>dt (66) 
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where p is a magnetic charge and the Dirac quantization condition is pe = 1. 
The usual skyrmion has a magnetic moment which would interact with a mag- 
netic monopole and break spherical symmetry. We use instead a non-topological 
ansatz for the chiral field 

U = e if{r ' t)r3 . (67) 

One can see that this field is made up of neutral pions and commutes with the charge 
matrix Q = \ + Nevertheless it has a non-zero total electric and baryonic charge 
due to the effects of anomalies as we shall see below. After inserting the ansatz (|67|) 
into the zeroth component of the baryonic current (|65|). one obtains 

n* = g[/(°o)-/(0)]. (68) 

The solution with the boundary conditions /(0) = and /(oo) = 2ir possesses 
unit baryon number and hence it can be interpreted as a baryon surrounding the 
monopole. 

If the field ansatz (J67|) is substituted into the meson and electromagnetic inter- 
action terms in the Lagrangian, they become 

C s = -^f{df? (69) 

C A = --^E i B i f (70) 

C em = ^(E 2 -B 2 ) (71) 

where the index i = 1,2,3 and the electromagnetic fields Ei and Bi are defined by 
Foi = y/—gtt Ei and Fij = CijkB k . When combined 

1 / p 2 \ 2 1 / p 4 



C A + £em _ _ - -Bf j - [1 -- 2 f<y (72) 

The extrema of the action occur when the electric field is given by 

E=^Bf. (73) 

This situation is reminiscent of the factorization of the Lagrangian that occurs for 
a BPS monopole |41| . 

The electric field implies a total charge 

or asymptotically 



Qoo = n B e 



(75) 
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and the ris = 1 solution can therefore be interpreted as a proton. If a black hole 
appears in the background, the inner boundary condition for the field / should be 
imposed not at the origin but at the event horizon r = r+. Thus the baryon number 
in the presence of an event horizon will be defined as 

n B = g[/(oo)-/(r + )], (76) 
which implies the baryon number swallowed by the black hole is 

n5 = g[/(r+)-/(0)]. (77) 

Thus the total baryon number can be recovered as the sum 

= n B + n B . (78) 

If /(r+) = 0, the baryon number is still an integer and conserved. This configuration 
represents a proton tightly bound to the black hole. On the other hand, if / takes 
some positive value at the horizon the baryon number is not an integer and the 
skyrmion carries fractional baryon number and electric charge. This configuration 
will be interpreted as a proton partially swallowed by the black hole. In particular, 
f(r+) = 27r means that the black hole has swallowed a whole proton, leaving nothing 
outside the horizon. 

It is interesting to observe that, while the baryon number disappears inside the 
horizon, the electric charge of the black hole can still be measured outside, turning 
the monopole black hole into a dyon black hole. Therefore, while the baryon number 
conservation is violated, charge conservation is not violated. 



4.1.1 Extremal black hole solutions 

In the extremal case we can obtain a general solution based on the Papapetrou- 
Majumdar metrics |421 14S| . We begin with the background metric fixed and later 
generalize to solve the full Einstein equations with chiral matter. The Papapetrou- 
Majumdar metrics have the form 

ds 2 = -H 2 dt 2 + H 2 {dx 2 + dy 2 + dz 2 ) (79) 

where 

H = l + ^GM 2L {m) 

n=l 

and R n is the ordinary Euclidean distance from the point mass M n located in three- 
dimensional space. We also associate these point masses with magnetic charges 
P n = G 1 ^ 2 M n , and the magnetic field 



(81) 
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The matter Lagrangian obtained earlier (|69|) has the form 



C = - l -Fl{dff-^B 2 (l + a 2 f) (82) 



where we will set 



a = ^, v 2 = ttGF 2 . (83) 

The Skyrme field equation obtained from the Lagrangian on this background be- 
comes 

- n 2 d 2 f + a 2 H- 2 (dH) 2 f = 0. (84) 

We lose no generality by taking equal charges P n = p. The solution with baryon 
number ng = n\j is then 

/ = 2ttH' s (85) 

where 



11a 2 

S = -2 + \/4 V (86) 

Since F n <C m p i, we can use s ps a/fj, for the pion model. 

For a single black hole, the Reissner-Nordstrom coordinate r is related to R by 
R = r — r + and we obtain 

Z = 2-(l-^f. (87) 

The field is effectively expelled by the black hole and vanishes on the horizon r = r + . 

The mass of the chiral field configuration can be obtained by integrating the 
Lagrangian (|%2*|) . 

m f = ^F 2 jjd l fdS* (88) 

where S is a large surface containing all of the masses. This gives 

m f = 27t 3 sF 2 GY J M n~K 3/2 n B eF w . (89) 

n 

The total mass in the chiral field is much less than one baryon mass per mass point. 
We can see how it is energetically favorable for a free skyrmion to change its internal 
configuration from the original Skyrme form to the simpler form used here when it 
comes into contact with a black hole monopole. The total topological description of 
this transformation for a single monopole is exactly as described in jUJ. 
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It is interesting to see that the electrostatic energy cancels due to the factor- 
ization occurring in the Lagrangian (|72|) . The chiral field mass is independent of 
the separation of the holes and therefore there are no forces between them. This 
is similar to the situation for BPS monopole solutions jlj, and suggests that there 
is a solution of the full Einstein-matter system. This existence of the solution will 
now be demonstrated. 

The spatial part of the Einstein tensor for the metric (|79|) is 

dj = -IH^im^djH) + H~ 2 (dH) 2 gij (90) 
and the Ricci scalar is 

R = -2H~ 1 d 2 H. (91) 

Substituting the Einstein tensor for the Lagrangian ()82|) into the Einstein equations 
gives 

H~ 1 d 2 H = -/i 2 (<9/) 2 (92) 
H~ 2 (diH)(djH) = -^(difXdjf) + GBiBjil + a 2 / 2 ) • (93) 

These make up a complete system of equations when we include the Maxwell equa- 
tion 



d(HB) = . (94) 

The second Einstein equation implies that df, dH and B are all parallel. We 
therefore impose a condition / = f(u), B = b{u)H~ 1 dH , where 

u = -^ 1 \ogH. (95) 

The system of equations becomes equivalent to an ordinary differential equation 
with independent parameter u, 

f" + Ml + f' 2 )f ~ u 2 b 2 f = (96) 

where 

,2 1 + f * 



1 + a 2 / 2 ' 



(97) 



The horizon corresponds to u — > — oo and the far region to u — > 0. The horizon 
must therefore be at a critical point of the first-order system corresponding to ()96|) . 
There is only one critical point, (/, /') = (0, 0), hence 

(/,/') -(0,0) as u^-oo. (98) 
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Since the critical point is a saddle, the solution is unique and exists for all values of 
\i. Having obtained the unique solution to ()96|) . we then define 

,o 

V(u) = l + n b(x)e~^dx . (99) 

J u 

It is easily seen from (|93)> that 

diV = V'diU = HB . (100) 
Hence the Maxwell equation implies d 2 V = and we can write 

V = 1 + ^^H. (101) 



n 



Inverting (|99jl gives u(V). 



4.1.2 Spherically symmetric solutions 

In the non-extremal case we shall consider spherically symmetric metrics which can 
be parameterized in the form 



d s 2 = -^e 25 dt 2 + '—dr 2 + rW (102) 
r z A 



where A and 5 are functions of r and t. After inserting the metric and the other 
field ansatz ()66|) and (|67|) into the Einstein field equations, one obtains 



\2 q^,4 ,>->4 ^-i4 

(Ae s f'y - -^e 5 f = --r~e- 5 Af - -e"^/ + -,e~ 5 } (103) 



8' = /A {^e~ 2S f 2 + / 2 ) (104) 

e-(^V = (105) 

where ^ and A are constants, 

^ = nF 2 G, A 2 = ^. (106) 

The electric charge within a sphere of radius r is given by equation IJ74JI . 

For very small p., which is the case for pions, the chiral field has little effect on 
the background metric and we may take 5 = and express A in terms of the mass 
M, electric charge Q and magnetic charge p of the black hole as 

A = r 2 - 2GMr + G{Q 2 + p 2 ) . (107) 
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The Skyrme field equation ()103j) on this background is therefore 

(A/')'-^/ = 0. (108) 

This should be solved subject to the boundary condition on which fixes the total 
charge, 

2 

Qoo = %r-foc = n B e . (109) 

The non-extremal black hole possesses two horizons at r = r_ and r = r + (r + > r_), 
related to the mass and charge by 

GM = ^(r- + r+) , GQ 2 = r^r + - Gp 2 . (110) 

The solution to Eq. (|108|) can be obtained analytically, 

P q {{r+ + r_)/(r+ - r_) - [2r + r_/(r + - r_)]/r) 
P s ({r+ + r_)/(r + - r_)J 

where -P s (2:) is a Legendre function and 



1 1 A 2 , . 
s = -- + W- + . 112 

2 y 4 r + r_ 

The black hole become a dyon with electric charge related to the value of / at the 
event horizon, 

n = V*l Mi3) 

W P s ((r + + r_)/(r + -r_))- { ) 

This relation can be solved, together with (|110j) . to obtain Q = Q(M), showing the 
existence of a one-parameter family of solutions (jib and p being regarded as fixed). 
In particular, Q — > as M approaches the extremal limit pM p i and the meson field 
is expelled from the hole. 

Larger values of [i may be realized for a hypothetical model where U is unrelated 
to pions, and this is discussed below. We will consider a static solution. We can 
replace A by a mass function m(r), defined implicitly by the relation 

A = r 2 - 2Gmr + G(p 2 + q 2 ) (114) 

where the charge is given by equation (|74|). The static equations become 

m ' = 2^/ + ^ A2// ' (U5) 

y = A(/') 2 (lie) 

f"+(^ + S')f'-^ K f = 0. (117) 
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Suitable boundary conditions are / — ^ 27r and 5 —* as r —* oo. 
In the numerical results the fields are scaled to the horizon size, 

f = — , rh = . (118) 

The solutions are parameterized by a parameter p, defined by 



P 2 = % (H9) 



which is restricted to p ^ 1. 

The extremal black hole solutions have A + = = 0. The regular solution to 
equation (|117|) has 

/+ = 0, Q = |-/ + = 0, p = l. (120) 

Hence the proton lies fully outside the black hole, as we saw before. The numerical 
solution for / is shown in Fig. E] This agrees with well with the result ()87|) , because 
the value of \x used here is still quite small. The results are still qualitatively similar 
for chiral models with \x of order one. 

For the non-extremal solution, we begin the integration of the field equations 
close to the horizon, with 

to = too + rh\{f — 1) + rri2{r — l) 2 + • • • (121) 
5 = 8 + + 8 x {f -!) + ■■■ (122) 
f = f + + f 1 (f-l) + ... (123) 

where 5+ and /+ are shooting parameters determined so as to satisfy the boundary 
conditions at infinity. 

As can be seen from the above expansion, the Skyrme field must have a nonzero 
value at the horizon otherwise the only allowed solution is the trivial one. Conse- 
quently the nonextremal black hole acquires an electric charge 

Q = ^U, (124) 

and allows the skyrmion to have fractional electric charge. The numerical results 
for this solution are shown in Figs. Eland El Again, these agree well with the fixed 
background for small values of [i. 

We have a single one parameter family of solutions with p < 1. In Figs. Eland 
El we plot the horizon radius r + and skyrmion mass mj as functions of black hole 
mass M. Figure El is related to the entropy of the black hole. The entropy of the 
black hole can be defined in the form 

S = ^A bh (125) 
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where Abh is the area of black holes, which is (vrr^.). The research in the entropy of 
black holes has revealed that the entropy of an extremal Reissner-Nordstrom black 
hole is zero, despite its finite size and it cannot be reached from the non-extremal 
state |44] . Thus the extremal black hole is a thermodynamically different object 
from the non-extremal one. 

The other figure shows how the proportion of the skyrmion which is swallowed 
by the black hole increases with the black hole mass. 

Figure shows how the horizon value of / changes as the coupling constant \x 
changes. As can be seen from the figure, for small (i the electromagnetic interaction 
is dominant so that the skyrmion is absorbed by the black hole to a lesser extent. 
On the other hand for large \x the gravitational interaction is dominant so that most 
of the skyrmion is absorbed by the black hole. 
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Figure 4: Profile function / as a function of r = r/r + for an extremal hole and 
H = 10" 4 . 

<5 




-2.5 J 



Figure 5: Backreaction S (xlO 3 ) as a function of r = r/r + for a non-extremal black 
hole, p = pm p i/r + = 0.9. Results for \i = 10^ 3 and n = are shown. 
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Figure 6: Mass function m as a function of r for p = 0.9. Results for fi = 10 3 and 
fx = 10~ 4 are shown 
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Figure 7: Horizon radius r + /M c as a function of the black hole mass M/M c for 
H = HT 4 . 




Figure 9: The value of / at the horizon for various values of fi. 
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4.2 Stability Analysis 



In this section we show that the skyrmion solutions which we have obtained are 
stable under spherically symmetric linear perturbations. We shall begin with the 
analysis of a skyrmion on the fixed background. 

In the fixed-background case the skyrme field is the only perturbed field and can 
be expanded about the skyrmion solution f by writing 



f(r,t)=f (r)+e^Z(r). 

Equation (|126|) is inserted into equation Q1U3|) with 5 
equation 



u 2 ^ 



(126) 

to obtain the eigenvalue 
(127) 



\2 „4 

-(AoO +Ak = -£- 

v ' r z A 
where the background equation has been used. 

If u is real and uj 2 > the solution is stable, and if u is imaginary and lo 2 < 
it is unstable since the mode can grow or decay exponentially under the small 
perturbation. To show which is the case we multiply both sides of equation (|127[) 
by £ and integrate in r from the horizon to infinity 

A 2 



^' 2 + 



dr 



Ar 



(128) 



where integration by parts and boundary conditions have been used. It can be seen 
that the integrands of both sides are positive definite, which means that uj 2 > 0. 
Hence the skyrmion on the fixed background is linearly stable . 

Next we analyze the stability of the skyrmion with backreaction. In this case we 
have to expand the metric as well as the skyrme field around the classical solutions 
fo, 5 and A 

f(r,t) = /o(r) + h(r,t) 
S(r,t) = 5 + <5i(r,t) 
A(r,t) = Ao + Ai(r,t). 

These are substituted into Eqs. (|1()^|) - (|1()5|) to obtain the following coupled equations 
up to first order 



Aofc/J + A f[ + Ai/£) e So 

5[ = 2/, 2 r/^/{ 
2/x 2 A 2 



A 2 



(<5i/o + /i)e' 



do 



An 



- S °fl 



-Mi + 



O 5o 



Equation (|131|) can be integrated with the help of the static field equation, 

Ax = -2/i 2 rAo^/i . 



(129) 
(130) 
(131) 

(132) 
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Substituting Eqs. (|13Uj) and (|132j) into Eq. H129j) one obtains the first order equation 
for fi 



2/i 2 (rA e 5 °^ 2 



+ 



A 2 (M C ) 2 5n 



h 



'A Q e 5 °f[ 

Setting /i(r, t) = £(r)e lw * one obtains an eigenvalue equation for £, 



5 



h 



2^(rA e^ 2 ) + 



A 2 (M C ) 



(133) 



(134) 



(135) 



We introduce the tortoise coordinate r* such that 

dr* _ 1 

Aoe 50 

and r* runs from — oo to +oo as r runs from r+ to +oo. Then Eq. (|134j) is reduced 
to the Sturm-Liouville equation 



dr 



*2 



where 



U 



2fi 2 



So 



A e s 



(136) 



(137) 



On the left-hand side we have r 4 , which makes the equation different from the 
previous eigenvalue equation with the fixed background. However, since r remains 
positive through the whole space, the same conditions for stability as the ordinary 
eigenvalue equation can be applied. As can be seen by examining U, U — > as 
r — > r + , (i.e. r* — > — oo), U — > C/qo > as r — > oo, and U > in between. 
In addition the solution does not change its shape for any value of the coupling 
constant /U. Therefore we can safely conclude that a skyrmion with backreaction is 
also linearly stable. 

When the solutions are stable, proton decay can only take place when extra 
particle fields are included in the model to carry the electric charge away. Since 
the underlying SU (5) theory does not require lepton charge conservation, it should 
satisfy AB = AL. The configuration we concern here is spherically symmetric and 
thus the theory is reduced to one space and one time dimension. In this case the 
bosonization technique can be applied in order to include leptons |45M46| . According 
the technique, we simply replace fermionic current to a real scalar field as 

1 



7T 



e^d v cj){r,t). 



(138) 



Then following the procedure described in Ref. we obtain the equations for the 
fixed background (approximately) 



(A/')' 



A 2 



(/-</>) 



A 

J a' 



f, 



(139) 
(140) 
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The stability arguments no longer apply. The dynamical process of a black hole 
swallowing a proton can be examined by solving these time-dependent field equations 
numerically. For a flat-space time, these were solved by Chemtob in Ref. [Uj where 
the catalysis cross section was estimated approximately to be lmb//3 with velocity 
(3, which confirms that the monopole catalysis proceeds at a strong-interaction scale 
without any suppression. 
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5 B = 2 Black Hole Skyrmions 

In this section black hole solutions of B = 2 skyrmions with axisymmetry are stud- 
ied. The study by Hartmann et al. showed the Einstein- Yang-Mills-Higgs theory 
possesses axially symmetric monopole and black hole solutions We follow 

their numerical technique to solve the Einstein-Skyrme model with axisymmetry. 
The obtained solution exhibits a torus shape in the energy and baryon density with 
a black hole at the center. Similarly to the case of B = 1, the baryon number is not 
conserved and takes fractional values outside the horizon. 

Recent studies of theories with large extra dimensions indicate that a true Planck 
scale is of order a TeV and the production rate of black holes massive than the Planck 
scale become quite large |50| I51j. Therefore, if a skyrmion represents a proton, this 
kind of black hole may be created in p-p collisions at the LHC in future. 

This work was carried out with the collaboration of T. Torii and K. Maeda. 



5.1 Axially Symmetric Configurations 

The Lagrangian for the Einstein-Skyrme system has been already given in Eq. IJ21JI . 
Since the B = 2 skyrmion is axially symmetric, we impose axially symmetric ansatz 
for the metric and chiral fields. We use the ansatz for the metric given in Ref. 



ds 2 



-fdt 2 + j{dr 2 + r 2 dd 2 ) + y 2 sin 2 9d(p 2 



(141) 



where / = f(r, 9), m = m(r, 9), and I = l(r, 9). 

The axially symmetric chiral fields can be parameterized by 



U = cos F(r, 9) + if ■ hr sin F(r, 9) 



(142) 



with % = (sin cos Tup, sin sin nip, cos 0) and = 0(r, 9). The integer n corre- 
sponds to the winding number of the Skyrme fields and for B = 2 we have n = 2. 
In terms of F and 0, the Lagrangian takes the form 



C s = cf+Cf 



(143) 



where 



C 



(i> 



C 



(2) 



8 m 

■2 



(d r F) 2 + ^(d e F) 2 + (d r G) 2 + ^(d e e) 2 sm 2 F 



+ 



n 



m 



r 2 s j n 2 | 



I 



2a 2 r 2 



id 



+ - 



n 2 m 



sm 



sin 2 sin 2 F 



(d [r Fd e] e) 2 + 



i 



n 2 m 



: (d e Fy} S m z e 



(d r e)2 + —(d e ey } sin 2 f sm 2 



sin 2 F 
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and we have defined the notation 

dtrFdmQ = d r Fd B Q - df>Fd r Q . 



(144) 



The baryon current in curved spacetime is defined in Eq. (|25[). The baryon number 
is then given by integrating B° over the hypersurface t = 0, 



B 



drdtidipy/ gW B° 

i J drdO (d [r .Fd e] 9)sin9(l - cos2F) 
- J dFdQ sin 9(1 - cos 2F) 



1 

2^ 



(2F-sin2F)cos9 



F u 0i 



F ,e 



(145) 



where (Fq, Bo) and (F±, Bi) are the values at the inner and outer boundary respec- 
tively. In flat spacetime we have 

(F 0) e ) = (tt,0) and (Fi.Gi) = (0,vr), 

which gives B = 2. In the presence of a black hole, the integration should be 
performed from the horizon to infinity, which change the values of Fq and allow the 
B to take fractional values of less than two. This situation can be interpreted as 
the black hole absorbing a skyrmion as was seen in the B = 1 case. 

The energy density is given by the zero-zero component of the stress-energy 
tensor 



8 m 



1 



1 



{d r FY + ^{d e FY + (d r Q) z + -^(deQY sin' F 



n 



r 2 sin 2 9 I 
n 2 m 



+ - 



8111 



sin F sin 



1 



+ 



1 f 
2a 2 r 2 



L 



(d [r Fd 0] ef 



n 2 m 



sin 



(d r F) 2 + ^(dgF) 2 } sin 2 Q 
(d r @) 2 + ^{d e @Y \ sin^Fsin 2 9 



sin 2 F . 



(146) 



5.2 Boundary Conditions 

Let us consider the boundary conditions for the chiral fields and metric functions. 
At the horizon r = r/j, the zero- zero component of the metric satisfies 



Sdo = -/(r h ,0) = O. 



(147) 



Noriko Shiiki and Nobuyuki Sawado 



33 



Regularity of the metric at the horizon requires 

m(r h ,9)=l(r h ,9) = 0. (148) 

The boundary conditions for F(r,9) and Q(r,9) at the horizon are obtained by 
expanding them at the horizon and inserting into the field equations which are 
derived from 5Cs/5F = and 5Cs/&® = respectively. Consequently we obtain 

d r F(r h ,9) = d r e(r h ,9) = 0. (149) 

The condition that the spacetime is asymptotically flat requires 

/ (oo, 9) = m (oo, 9) = I (oo, 9) = 1 . (150) 

The boundary conditions for F and at infinity remain the same as in flat spacetime 

F(oo,0) = O, d r @{oo,9) = 0. (151) 

For the solution to be axially symmetric, we have 

(r, 0) = d e m (r, 0) = dgl (r, 0) = 0, (152) 
r, J) = d e m (r, J) = d e l (r, J) = 0. (153) 



, 2 j -v».y, 2 

Regularity on the axis and axisymmetry impose the boundary conditions on F and 
as 

d e F(r,0) =d g F(r,^j = 0, (154) 
9(r,0)=0, G(r,|) = |. (155) 

Under these boundary conditions, we shall solve the Einstein equations and the 
matter field equations numerically. 



5.3 Numerical Results 

Let us introduce dimensionless coordinate and coupling constant 

x = aF n r , a = TrGf% . 

Then the free parameters are the horizon Xh and the coupling constant a for this 
system. We shall take a = as decoupling of gravity from the matter, effectively 
G = 0. 

In Figs. ITUl 1111 are the energy densities of the black hole solutions with a = 0.0, 
1.5 respectively. As a becomes larger, the energy density becomes smaller and 
sparse. This can be interpreted that the black hole absorbs more skyrmions for a 
larger coupling constant. The shape is slightly distorted in the background of the 
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Figure 10: The energy density e in cylindrical coordinates p and z with Xh = 
1.0, a = 0.0. 

black hole so that one can see the spherically symmetric horizon at the center of the 
skyrmion. Fig. ^] is the baryon density around the black hole. The dependence of 
the baryon density on the value of the coupling constant is small. It can be checked 
that the energy and baryon density vanish at p = 0. Inserting the metric functions 
as well as the profile functions expanded around the horizon instead, one can also 
see that the energy and baryon density vanish at the horizon. 

The domain of existence of the black hole solution is shown in Fig. 1131 There 
exist minimum and maximum value of and a beyond which no black hole solutions 
exist. Therefore the regular skyrmion solutions can not be recovered from the black 
hole solutions by taking the limit of Xh — ► unlike the case of B = 1 |17j . In Fig. ll4l is 
the dependence of the baryon number on x^. One can see that the baryon number 
decreases as the black hole grows in size. This figure confirms that the baryon 
number is no longer conserved due to the black hole absorbing the skyrmion. 
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Figure 11: The energy density e in cylindrical coordinates p and z with Xh = 
1.0, a = 1.5. 




Figure 12: The baryon density 6 with xt = 1.0, a = 0.0. 
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Figure 13: The domain of existence of the solution. For a > 2.0, there exists no 
non-trivial solution. 




Figure 14: The dependence of the baryon number on the size of the horizon. 
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6 Conclusion and Discussion 

In this paper black hole solutions with Skyrme hair are reviewed. Sec. 3 is devoted 
to the study of the spherically symmetric black hole with B = 1 Skyrme hair and 
its stability analysis. In Sec. 4 the gauged Einstein-Skyrme system is constructed 
and the monopole black hole with B = 1 Skyrme hair is obtained. The extended 
Einstein-Skyrme system to B = 2 axially symmetric configuration is studied in 
Sec. 5. The black hole solution with B = 2 Skyrme hair is obtained and it is 
shown that the energy and baryon density exhibit a torus shape with the spherically 
symmetric horizon at the center. 

The common feature in those solutions is that they can support fractional bary- 
onic charges outside the horizon, violating baryon number conservation. The study 
of the monopole black hole skyrmion shows that although global charge conserva- 
tion such as baryon number is violated, gauge charge conservation such as electric 
or magnetic charge is still hold. It is remarkable that the black hole and monopole 
black hole solution with B = 1 Skyrme hair turned out to be stable under linear 
perturbations even with the non-integer baryonic charge. Obviously it is important 
to study the stability of the B = 2 solutions. We expect that the stability analysis 
may be performed by applying the catastrophe theory for black holes with non-linear 
hair [52]. 

For these microscopic black holes, however, we cannot ignore the quantum ef- 
fects. In fact, they are stable only classically and will decay to the unstable solutions 
due to quantum transitions. Besides black holes of the size of a proton should have 
large fluxes of Hawking radiation j^Uj ■ Hence situations in which baryon decay pro- 
cess might become more realistic and significant over the Hawking radiation occurs 
only when the black hole carries electric or/and magnetic charge with which the 
skyrmion interacts electromagnetically as well as gravitationally. Especially inter- 
esting is the extremal black hole which has a vanishing effective temperature, so 
the Hawking radiation may even vanish. The free skyrmion has a magnetic moment 
and, if it has the correct orientation, it will be attracted to the monopole. When the 
proton approaches the black hole monopole, the fields rearrange themselves into the 
energetically preferred configuration of skyrmion hair solutions described in Sec. 4. 
For the stable solutions, baryon decay can only takes place when extra particle fields 
are included. We have given a rigorous argument to study dynamical proton decay 
to a lepton by a monopole black hole. The bosonization technique is conveniently 
introduced and the system is reduced to two coupled time-dependent field equations. 
Solving these equations with appropriate boundary conditions at the horizon will 
be our future work. 

Recently we performed collective quantization of a B = 1 gravitating skyrmion 
and calculated various nucleon observables under the gravitational effects in Ref. |53| . 
It will be interesting to extend the work to the black hole skyrmion. Also interest- 
ing but hard is to consider black holes with B > 3 Skyrme hair which has discrete 
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symmetries. The rational map ansatz will be a powerful tool for obtaining multi- 
skyrmion black hole solutions. The inclusion of gauge fields may also be possible to 
study the interaction between a monopole black hole and multi-skyrmion. Extend- 
ing the model to higher dimensions will be also an exciting problem. It is known 
that multi-extremal black hole solutions can be generalized to p-brane objects by 
coupling to an antisymmetric tensor A p +i with p + 1 indices. In brane theory con- 
text, Skyrme fields are axions interacting with branes. The generalization of our 
model to brane theory will be worth studying in future. 
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